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Abstract _ _ _ 

The  front-tracking  “hydrocode”  FRONTIER  limits  finite-difference  solution  of  the  dynamic 
continuum  mechanics  equations  to  regions  bounded  by  tracked  fronts  (shock  waves,  contact 
discontinuities,  and  gradient  discontinuities)  and  uses  local  solutions  of  the  Riemann  problem 
to  advance  the  positions  of  these  fronts.  This  solution  method  places  stringent  requirements  on 
the  availability  of  thermodynamics  information.  With  the  addition  of  an  appropriate  reacting 
mixture  model,  the  code  can  be  applied  to  solid  explosive  initiation  problems.  The  required 
thermodynamic  functions  are  more  difficult  to  obtain  for  reacting  mixtures  even  though  the 
equations  of  state  for  each  of  the  phases  present  are  known.  W e  have  developed  a  mixture  model 
based  on  assumptions  of  mechanical  equilibrium  and  thermal  isolation  that  can  be  used 
independent  of  the  choice  of  equations  of  state  for  the  phases,  and  we  have  derived  expressions 
for  the  necessary  thermodynamic  functions. 
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1.  Background 


We  are  working  to  adapt  FRONTIER,  a  front-tracking  "hydrocode"  developed  at  the  State 
University  of  New  York  at  Stony  Brook  (Glimm  et  al.  1981,  1985;  Grove  1989,  1993),  for 
application  to  solid  explosive  initiation  problems.  The  principal  elements  required  to  treat  such 
problems  numerically  are  a  continuum  mechanics  model  (including  appropriate  treatment  of 
boundary  conditions)  and  a  reacting  mixture  model.  The  former  is  represented  by  a  system  of 
dynamic  partial  differential  equations  and  the  latter  by  a  system  of  algebraic  equations  (augmented, 
in  some  cases,  by  additional  dynamic  equations).  The  forms  of  these  models  depend  on  the 
assumptions  used. 

Frontier  is  unique  in  that  it  limits  finite-difference  solution  of  the  dynamic  continuum 
mechanics  equations  to  regions  bounded  by  tracked  fronts.  These  may  include  shock  waves,  contact 
discontinuities  (i.e.,  material  interfaces  and  sliplines),  and  gradient  discontinuities  (e.g.,  boundaries 
between  simple  and  complex  wave  regions).  The  front-tracking  method  superimposes  grids  for  these 
fronts  on  the  computational  region.  For  an  n -dimensional  computational  region,  the  grids  are  of 
dimension  n-1.  The  method  uses  local  solutions  of  the  Riemann  problem  coupled  to  the  interior 
region  solution  to  advance  the  positions  of  tracked  shock  waves  and  contact  discontinuities.  There 
results  a  description  of  the  flow  that  retains  the  discontinuous  character  of  the  fronts.  This  solution 
method  places  stringent  requirements  on  the  availability  of  thermodynamic  information.  In  addition 
to  the  incomplete  equation  of  state  for  pressure  required  by  all  hydrocodes,  the  Riemann  solvers 
require  expressions  for  the  acoustic  impedance,  sound  speed  and  adiabatic  exponent,  and  the 
Griineisen  function,  as  well  as  adiabats  and  Hugoniots.  While  these  are  readily  developed  from  most 
equations  of  state,  they  are  more  difficult  to  obtain  for  reacting  mixtures. 

In  order  to  avoid  the  need  to  completely  rework  FRONTIER'S  coding  for  solution  of  material 
dynamics,  we  are  interested  in  mixture  models  that,  insofar  as  possible,  allow  retention  of  the  gas 
dynamics  conservation  equations  for  diffusionless  flow.  The  extent  to  which  this  objective  may  be 
achieved  depends  on  the  physical  assumptions  that  we  wish  to  apply.  Usually,  the  requirement  can 
be  met,  although  additional  dynamic  equations  may  arise.  It  is  generally  necessary  to  add  a  dynamic 
equation  representing  reaction  progress. 
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Several  modeling  approaches  applicable  to  mixtures  of  discrete  phases  have  been  reported  in  the 
literature,  and  some  of  these  are  reviewed  and  amplified  here.  The  set  of  physical  assumptions  that 
appears  to  best  apply  to  explosive  initiation  includes  mechanical  equilibrium  and  thermal  isolation 
between  reactant  and  product  phases.  While  the  thermal  isolation  condition  can  be  implemented  as 
an  added  dynamic  equation,  representing  partitioning  of  the  energy  conservation  equation  between 
the  phases,  including  it  in  the  mixture  model  is  more  consistent  with  our  desire  to  retain  the  original 
form  of  the  conservation  equations.  We  have  developed  a  model  along  these  lines  that  can  be  used 
independently  of  the  choice  of  equations  of  state  for  the  phases  and  have  derived  expressions  for  the 
necessary  thermodynamic  functions. 

2.  Mixture  and  Phase  Variables 

Consider  a  Lagrangian  control  volume,  A  V,  containing  a  mixture  of  two  phases  having  mass  AM 
Phase  variables  are  denoted  by  subscripts  r  and  p  (anticipating  but  not  limited  to  reactants  and 
products),  while  mixture  variables  are  unsubscripted.  The  mass  of  the  mixture  is  equal  to  the  sum 
of  the  phase  masses: 

AM  =  AM^+ AMp. 

Dividing  by  AM  gives 

AM,  ,  AM^ 

AM  AM 

Mass  fractions  for  each  phase  can  be  defined  as 

AM, 

A,=  lim  - , 

AM-O  AM 
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and 


A  =  lim  - ^  , 

A/w-o  AM 


so  that 


A^+ Ap-  1 


In  light  of  this  simple  relation,  a  single  mass  fraction  variable,  A,  may  be  chosen  so  that 


A,=  1-A 


Ap=A, 


Similarly,  volume  fractions  are  defined  as 


nr=1-n=  lim  —  , 
AV-o  AV 


np=n=  iim 

^  AV-o  AV 


The  mixture  and  phase  specific  volumes  are  defined  as 


AV 
lim  - 

AM-O  AM 
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LVr 

v=  iim  — — 

LM-O  AM^ 


and 


A\/„ 


v=  Iim  — 

^  AM-O  AM„ 


Thus, 


(1 


and 


KVp=r\v. 


The  mixture  is  said  to  be  saturated  if  its  volume  equals  the  sum  of  the  phase  volumes. 

AV  =  AV^+AVp. 

This  assumption  is  most  appropriate  when  one  of  the  phases  is  gaseous,  as  is  generally  the  case  for 
detonation.  Dividing  by  AM  gives 

^  _  AV^  ^  ^  AVp  AMp 
AM  AM,  AM  AMp  AM 

so  that 


Similar  considerations  yield  samration  conditions  for  the  mixture  and  phase  specific  internal  energies 
and  specific  entropies: 
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e  =  (1  -A)e,+  Aep, 


and 


s  =  (1  -K)s^+KSp. 


3.  Conservation  Equations 

It  is  assumed  that  the  values  of  some  of  the  variables  required  to  describe  the  mixture  are 
available  from  integration  of  continuum  mechanics  conservation  equations.  Since  the  material 
strengths  of  explosives  are  generally  modest  and  the  pressures  associated  with  their  initiation  high, 
the  equations  of  gas  dynamics  can  be  used: 

+V'(p£/)  =0, 

dt 

—  (pu)  +V-(pi/i/)  +Vp  =  0, 


and 

^{pej)  +  V-ipejU)  +  V-(pi/)  =  0, 

dt 

where  t  is  time,  p  is  the  density,  u  is  the  particle  velocity  vector,  p  is  the  pressure,  e-j-  is  the  total 
specific  energy, 

er=  e*-U’U, 

'  2 
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and  e  is  the  specific  internal  energy.  This  system  is  generally  closed  by  obtaining  pressure  as  a 
function  of  density  (or  specific  volume,  v'-  1/p)  and  specific  internal  energy  from  an  incomplete 
equation  of  state: 

p  =  p{v,e). 


In  order  to  apply  the  gas  dynamics  equations  to  a  discrete  mixture,  an  incomplete  equation  of 
state  yielding  a  single  pressure  for  the  mixture  is  required.  This  pressure  will,  additionally,  depend 
on  a  reaction  progress  variable: 

p  =  p{v,e,K). 

Further,  the  phase  pressures  should  equal  the  mixture  pressure.  Otherwise,  the  momentum  and 
energy  conservation  equations  must  be  split  between  the  phases.  Splitting  is  also  required  if  the  flow 
is  not  diffusionless  (i.e.,  the  same  particle  velocity  should  apply  to  both  phases). 


Where  the  mixture  is  reacting,  the  reaction  progress  must  be  computed,  and  the  energy  liberated 
during  the  reaction  must  be  accounted  for.  Additional  variables  and  equations  are  required  for  these 
purposes.  However,  the  conservation  equations  in  the  form  shown  apply  as  long  as  single  pressures 
and  particle  velocities  are  maintained.  The  reaction  rate  equation  is 

—  (pA)  +  V-{pAu)  =  pA(A,p . ) , 

dt 

where  the  reaction  rate  is  a  known  function  of  state  and  history  variables  obtained  from  one  of  a 
number  of  available  models.  The  equation  arises  as  a  consequence  of  conservation  of  mass  and  must 
be  integrated  along  with  the  other  conservation  equations. 

The  chemical  energy  released  in  the  reaction  process  must  also  be  accounted  for.  This  could  be 
done  by  modifying  the  energy  conservation  equation  such  that 
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—  (pe7-)  +  V'(pe7-u)  +V-(pu)  =  -pAh[jk{K,p . ), 


where  is  the  (negative)  heat  of  detonation.  However,  this  is  inconsistent  with  our  desire  to 
retain  the  original  conservation  equations  and  it  can  be  accomplished  more  simply  by  regarding  the 
reactant  phase  specific  internal  energy  as  composed  of  a  thermal  component  and  a  potential 
component  equal  to  the  absolute  value  of  the  heat  of  detonation: 

~  thermal  ~  ' 

The  reactant  phase  equation  of  state  (as  developed  in  Appendix  A)  applies  only  to  the  thermal 
component, 

^r,  thermal  ~  ^r'^  ’ 

and  the  reactant  phase  equation  of  state  for  pressure  may  be  expressed 

p^=  p^{v^,e+LhQ) . 

In  practice,  the  heat  of  detonation  is  often  included  in  the  equation  of  state. 

4.  Mixture  Equations  of  State 

4.1  General.  Returning  to  our  Lagrangian  control  volume,  we  may  write  the  complete  equation 
of  state  for  the  mixture: 

AE=E(Al/,AS,AM,,AMp), 

where  AS  is  the  entropy.  Because  of  the  homogeneous  first-order  nature  of  such  equations,  dividing 
by  AM  gives 

AE.  =  e{—  —  ^  ' 

AM  AM  ’  AM  ’  AM  ’  AM 
\  / 
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Thus,  in  the  limit  as  AM  vanishes, 


Q  —  S,A^,  Ap  )  —  6{VfS,\) 


Incomplete  equations  of  state  are  obtained  by  differentiating: 

p{v,S,K)  =  -{ 


'de^ 


dv 


s.X 


and 


r(i',s,x)  = 


f 


where  T  is  the  temperature.  Additionally, 

/  V \  f  de^ 

V  dA  / 

defines  the  normalized  chemical  potential.  By  solving  the  complete  equation  of  state  for  the  specific 
entropy, 

S  =  s(</,e,A) , 


and  substituting  into  the  incomplete  equation  of  state  for  pressure,  we  obtain  the  form  required  for 
solution  of  the  gas  dynamics  conservation  equations  (i.e.,  with  v,  e,  and  A  as  the  three  independent 
variables): 

p  =  p{v,e,K) . 


However,  the  complete  equation  of  state  is  generally  unavailable.  Rather,  incomplete  equations 
of  state  for  each  phase  giving  pressure  as  a  function  of  specific  volume  and  internal  energy  are 
known: 

p^=  p^(v;.,e^+Afi£,) , 


and 
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Pp  =  Pp{Vp,ep). 

In  order  to  establish  an  equation  of  state  for  the  mixture,  we  may  combine  these  with  the  volume  and 
energy  saturation  conditions: 

v={^  -\)v^+\Vp 

and 

e  =  (1  -A)e,+  Aep. 

This  is  a  system  of  four  equations  in  nine  variables  (v,  e,  A,  ,  p^ ,  Vp  ,  Bp  ,  and  Pp ),  not 

including  the  pressure.  Thus,  two  additional  equations  are  required  to  produce  a  solution  for  the 
pressure  as  a  function  of  three  independent  variables. 

4.2  Equilibrium  Assumptions.  The  number  of  independent  variables  can  be  reduced  by 
invoking  assumptions  of  physical  equilibrium.  The  most  common  of  these,  the  assumption  of 
mechanical  equilibrium  (Mader  1979;  Johnson,  Tang,  and  Forest  1985;  Starkenberg  1989),  is 
required  to  retain  the  simplified  form  of  the  conservation  equations.  It  implies  that  the  phases  are 
at  equal  pressures  and  is  achieved  by  a  relatively  rapid  wave  propagation  mechanism.  Thus, 

p^{v^,B+Ahu)=Pp{Vp,Bp). 


Another  assumption  that  has  been  used  (Mader  1 979)  is  that  of  thermal  equilibrium.  In  this  case, 
the  phases  are  at  equal  temperatures,  presumably  by  virtue  of  significantly  slower  conductive 
processes,  and  incomplete  equations  of  state  for  the  phase  temperatures  are  required: 

Tr{v„B+Ahij)  =  Tp{Vp,Bp). 
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Thus,  with  both  mechanical  and  thermal  equilibrium,  the  mixture  equation  of  state  may  be 
expressed  as 

Vp{v,e,\), 

e=e^{v,e,K), 

ep=ep{v,e,K), 

p  =  p{v,e,\), 

T=T{v,e,k). 

This  formulation  has  the  desired  set  of  independent  variables.  However,  because  thermal 
equilibrium  can  be  achieved  only  through  physical  processes  that  generally  act  slowly  compared  to 
pressure  equilibration  and  are  consequently  physically  unreaUstic  in  explosive  initiation 
environments,  alternatives  to  this  assumption  have  been  sought. 

4.3  Thermal  Isolation.  The  thermal  isolation  condition,  which  represents  the  opposite  extreme 
and  appears  more  realistic  than  thermal  equilibrium,  was  used  in  Lagrangian  computations  by 
Johnson,  Tang,  and  Forest  (1985)  and  in  conjunction  with  smoothed  particle  hydrodynamics  by 
Libersky  (1993).  These  approaches  were  limited  to  use  with  a  linear- f  Mie-Gruneisen  equation  of 
state  for  the  reactant  phase.  Here,  we  shall  relax  this  limitation  and  generalize  the  approach  for  use 
with  any  equation  of  state. 

If  the  mixture  is  constrained  to  change  its  state  such  that  the  phases  remain  thermally  isolated, 
there  results  an  isentropic  process  for  the  reactant  phase.  Thus,  the  thermal  isolation  condition  may 
be  expressed 

Pr{V^,e+Ah[y)  =  PsMr’'^rk^^rk^^^D)  > 
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or 


e,^^h^  =  esr{v,;v^,e^^Mi^), 

or  simply 

SMr,  V  A/7d)  =  Sr('^rfc-  > 

where  the  isentrope  functions  are  known  and  (v^f,  e^^f)  is  a  reference  state  specifying  the  isentrope. 
If  the  material  has  been  shocked,  this  state  is,  most  conveniently,  that  existing  on  the  Hugoniot 
immediately  after  the  most  recent  shock  passage.  Otherwise,  the  initial  state  may  be  used.  In  an 
Eulerian  numerical  computation,  this  state  can  be  advected  for  use  throughout  the  reactive  flow  field: 

Ol 

and 

|-(pej  +  v-(pe*u)=0. 

at 

Unfortunately,  where  artificial  viscosity  is  used  to  capture  shock  waves,  the  accurate  determination 
of  shocked  states  (and,  hence,  the  accurate  initialization  of  the  reference  state  following  shock 
passage)  is  not  generally  possible.  The  front-tracking  approach,  on  the  other  hand,  provides  an 
opportunity  to  accurately  initialize  and  advect  Srfc)- 

4.4  Other  Models.  Because  the  iterative  processes  associated  with  the  foregoing  approaches 
may  be  time  consuming  (although  convergence  should  generally  be  rapid),  an  alternative  approach 
(Kerley  1992)  has  been  used  in  the  Sandia  National  Laboratory’s  CTH  code.  This  is  referred  to  as 
the  "two-state"  model.  Here,  the  mixture  equation  of  state  is  expressed  as  follows: 

p{p,r,A)  =  {i-A)p,(p,r)+App(p,r)  , 


and 
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e(p,7-.A)  =  (1  -A)  e/p.T)  +  A  ©^(p.r)  • 


Thus,  each  phase  is  assumed  to  have  the  density  and  temperature  of  the  mixture,  and  the  mixture 
pressure  is  a  weighted  average  of  the  phase  pressures  and  has  no  physical  meaning.  Because  of  the 
poor  physical  fidelity  associated  with  these  assumptions,  we  did  not  consider  it  to  be  a  desirable 
candidate  for  incorporation  into  FRONTIER.  The  approach  has  been  eliminated  for  the  1999  release 
ofCTH. 

In  addition,  Baer  and  Nunziato  (1986)  have  developed  a  nonequilibrium  model  with  diffusion 
for  CTH  based  on  physical  assumptions.  Because  this  model  was  primarily  intended  for  porous 
materials,  has  exhibited  difficulties  (Menikoff  1996),  and  requires  substantially  more  complex 
conservation  equations,  it  was  not  considered  for  application  in  the  front-tracking  environment. 

4.5  Physical  Considerations.  To  this  point,  our  physical  assumptions  have  been  driven  by  the 
heuristic  purpose  of  retaining  the  gas  dynamics  conservation  equations.  In  spite  of  this,  there 
remains  a  variety  of  circumstances  under  which  the  assumptions  of  diffiisionless  flow  in  mechanical 
equilibrium  and  thermal  isolation  may  be  regarded  as  physically  reasonable. 

The  physical  justification  for  these  assumptions  begins  with  the  notion  that  the  process  of 
explosive  initiation  starts  at  hot  spots  that  are  excited  by  the  interaction  between  a  passing  shock 
wave  and  the  heterogeneities  in  the  explosive.  Such  hot  spots  have  been  observed  experimentally 
by  von  HoUe  and  Tarver  (1981).  They  are  distributed  on  the  scale  of  the  heterogeneities  and 
represent  a  very  small  fraction  of  the  explosive  mass.  There  are  assumed  to  be  numerous  hot  spots 
in  any  control  volume  used  to  derive  equations  of  motion.  The  hot  spots,  which  may  be  idealized 
as  closed  regions,  are  assumed  to  react  thermally  and  to  produce  "holes"  bounded  by  closed  surfaces. 
The  reaction  is  assumed  to  propagate  further  by  a  surface  burning  mechanism  akin  to  that  observed 
for  propellants  in  strand  burner  experiments,  which  occurs  at  uniform  pressure.  This  portion  of  the 
reaction  process  is  often  referred  to  as  hole  burning.  As  reaction  proceeds,  the  holes  join,  and,  before 
long,  discrete  solid  particles  appear  embedded  in  a  gas  phase  matrix.  At  this  point,  the  process  is 
called  grain  burning.  The  point  of  transition  from  hole  burning  to  grain  burning  depends  on  the 
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shape  of  the  holes  and  the  way  in  which  they  evolve.  This  description  is  applicable  to  the  initiation 
of  condensed  energetic  material  by  a  process  referred  to  as  shock-to-detonation  transition.  Porous 
materials  exhibit  a  more  complex  initiation  process  referred  to  as  deflagration-to-detonation 
transition  in  which  diffusion  may  play  a  role. 

The  assumption  of  diffusionless  flow  is  easiest  to  defend  within  the  surface-burning  context. 
During  the  early  stages  of  reaction,  the  gas  is  entrapped  within  the  solid  phase  and  may  not  diffuse. 
During  the  late  stages  of  reaction,  high  drag  forces  acting  on  the  remaining  solid  resist  any  tendency 
toward  diffusion.  But,  most  importantly,  since  burning  is  assumed  to  occur  at  closed  surfaces 
throughout  the  process,  the  net  momentum  of  the  products  relative  to  the  reactants  tends  to  vanish. 
Mechanical  equilibrium  is  more  difficult  to  justify.  It  rests  on  the  assumption  that  burning  occurs 
as  observed  on  a  macroscopic  scale  for  propellants.  Assuming  thermal  equilibrium  appears 
unjustified  for  the  rapid  processes  associated  with  initiation  of  detonation.  Thermal  isolation  lies 
at  the  opposite  extreme  and  is  valid  when  the  reaction  process  is  sufficiently  rapid. 

Wackerle  and  Anderson  (1984)  have  shown  that  simplified  representations  of  the  surface  reaction 
rate  and  the  equilibrium  state  may  be  inconsistent  with  the  assumed  reaction  topology.  However, 
in  the  context  of  an  explosive  initiation  model,  all  of  this  may  be  moot.  Such  a  model  is  completed 
by  specifying  a  reaction  rate  function  containing  constants  that  are  calibrated  with  respect  to  certain 
explosive  sensitivity  experiments.  Thus,  some  excellent  results  are  achieved  with  highly 
questionable  "equilibrium"  assumptions.  No  advantage  to  more  complex  models  has  yet  been 
exhibited. 

5.  A  Discrete  Mixture  Equation  of  State 


5.1  General.  The  equations  governing  the  mixture  in  mechanical  equilibrium  and  thermal 
isolation  are 
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i^=(1  -K)v^+KVp, 

e  =  (1  -\)e^+Kep, 

P= 

P=Pp('^p'®p)' 

and 

P=  PsriV>^rk’eri<^^^D)’ 

or 

e/A/7o  =  65,(1/,;  V,e^+A/7o). 

This  constitutes  an  incomplete  equation  of  state  for  the  mixture  pressure  and  phase  states  that 
includes  the  reactant  reference  state  as  an  independent  variable: 

v^=  v^{v,e,K]v^,e^+Ah[)), 

^p=  Vp{v,e,K’,v^,e^+Lhc), 
e^=  eAv,e,K;v^,e^+Ah[j) , 
ep=ep{v,e,K-,v^,e^+Ah^), 
and 

p  =  p{v,  e,\;  e^+Lhg) . 


5.2  Mixture-Equation  Solutions.  Because  closed  form  solutions  to  the  system  of  equations 
describing  the  mixture  cannot  generally  be  found,  it  is  necessary  to  obtain  solutions  using  iterative 
approaches.  It  is  assumed  that  initial  estimates  for  the  values  of  the  unknown  variables  are  available: 


,(0)  „(0)  -(0)  -(0) 
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If  these  are  known  to  sufficient  accuracy,  successively  more  accurate  solutions,  given  by 


and 

may  be  obtained  using  Newton's  method.  Other  methods  may  be  used  until  accuracy  sufficient  for 
convergence  of  Newton's  method  has  been  achieved.  Substituting  these  solutions  into  the  original 
equations  (eliminating  the  superscripts  for  clarity)  yields  the  variational  equations. 


(1  -K)dv^+K5Vp=  i^-(1  -K)v^-KVp, 

(1  -A)56j.+  A66p=  6  ~  (1  “A)  6^.-  ASp, 


bv  + 

'Sp,' 

5©  - 

5v.- 

- ^ 

CD 

r 

Or 

r 

Vr 

P 

en 

CD 

6ep  =  Pp(Vp , Bp)  - e,+ Aho) , 


and 


PsriV’Vrto^rk^^l^D)  5©,  =  CsriV/,  ,  6^^+ A/Iq)  -  (e/A/)o)  . 


These  equations  are  then  solved  for  ,  bVp  ,  6©^  and  5  ©p,  and  the  solution  is  updated.  The 
procedure  is  repeated  until  the  desired  accuracy  is  achieved. 
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Li  hydrocode  applications,  and  can  be  evaluated  by  calls  to  appropriate  equation-of-state 
subroutines.  The  isentrope  reference  state,  on  the  other  hand,  is  an  integrated  variable  of  the  flow 
and  must  be  included  in  the  call  to  the  mixture-equation-of-state  routine. 


5.3  Mixture  Derivatives.  Expressions  for  the  derivatives  of  the  independent  state  variables 
with  respect  to  the  dependent  state  variables  are  developed  in  Appendix  B.  Derivatives  of  the 
pressure  are 


_  1 

and 


/  3., 

dp 

if 

^Pp\ 

<  dej 

where 

and  //  and  ip  are  the  squares  of  the  acoustic  impedances  of  the  phases,  defined  as 


dp 


m 


dV. 


=  P 


I  ^Pm 


m  ) 


m 


de 


m 


dv, 


m  =  r,p 


m  ) 


With  these  derivatives  in  hand,  we  may  proceed  to  relate  thermodynamic  functions  of  the  mixture 
to  their  phase  counterparts. 

5.4  Sound  Speed  and  Related  Functions.  The  square  of  the  acoustic  impedance  of  the  mixture 
may  be  written 
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i2  =  JiP 


dvL,  laej„x  lav^Jpx 


It  varies  continuously  from  that  of  the  reactants  to  that  of  the  products  as  the  mass  fraction  varies 
from  0  to  1,  and  its  reciprocal  is  subject  to  a  saturation  condition: 


J_  =  +A 


The  square  of  the  sound  speed  may  be  defined  in  terms  of  the  square  of  the  acoustic  impedance: 


{1-K)v^c^^hv^cf 


The  corresponding  saturation  condition  is 


Another  important  function  defined  in  terms  of  the  square  of  the  acoustic  impedance  is  the 
adiabatic  exponent,  given  by 

=  ^VrYp  ^  ^r^p 

P  (1-A)v^,Yp+AV^  (1-n)Yp+nYr 

with  the  saturation  condition 


X.  =  (1-A)^  +  A-^ 
Y  Y,  Yp 
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5.5  Griineisen  Function.  The  Gruneisen  function  of  the  mixture  is 


.I^r-Lr 

"  Y,  ' 


The  continuity  of  this  expression  for  vanishing  X  has  not  been  established. 

5.6  Temperature  and  Specific  Heat.  In  order  to  determine  the  thermodynamic  temperature 
of  the  mixture,  we  begin  by  differentiating  the  saturation  condition  for  the  entropy  with  respect  to 
S: 


I 


+  K 


/ 


dS 


=  1 


Noting  that  the  thermal  isolation  condition  implies  that  is  a  constant,  we  find  that 

(dsA 


ds 


=  0, 


and 


(  dSp] 

ds 


v;A 


X 


Differentiating  the  saturation  condition  for  the  specific  volume  with  respect  to  S  gives 


I 


+  A 


v,X 


( 

ds 


=  0. 


v,X 


We  may  write  the  mechanical  equilibrium  condition  in  the  form 

p,{v^,s;^=pp{Vp,Sp), 

and  differentiate  with  respect  to  S  to  obtain 
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dS 


-  I 


v,k 


v,X 


-JVfi 


Thus, 


and 


dv,] 


v,A 


[pIe 


f 


\ 


\ 


) 


v,A 


^[pjp  ^ 

K  v^P 


Finally,  the  temperature  is  determined  from  its  definition, 


T  = 


dS 


by  differentiating  the  saturation  condition  for  the  specific  internal  energy. 


e{v,s,K)  =  0-K)  e,(v^  ,s,)  +  K  ep{Vp ,  s^) . 


This  shows  that  the  mixture  temperature  is  simply  equal  to  the  product  temperature: 


T{v,e,K)  =  Tp[Vp(v,e,K),ep{v,e,K)]. 
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The  latter  can  be  computed  if  an  incomplete  equation  of  state  for  the  products  giving  temperature 
as  a  function  of  specific  volume  and  specific  internal  energy  is  available.  The  continuity  of  this 
expression  for  vanishing  A  has  not  been  established. 

Similarly,  the  specific  heat  at  constant  volume  of  the  mixture  may  be  determined  from  its 
definition, 


by  differentiating  the  expressions 


e{v,  T,h)  =  (1  -A)  e,{v^,  7,)  +  Aep(Vp ,  T^) , 

v=0  -K)v^+KVp, 


and 


=  PsMr)  ’ 


with  respect  to  T.  This  shows  that  the  mixture  specific  heat  is  given  by 


A /2  c., 


(  f 

/2-(i-A)c  r„  -£ 

"  \''pi 


Similarly,  the  continuity  of  this  expression  for  vanishing  A  has  not  been  established. 


5.7  Adiabat  The  complete  equation  of  state  of  the  mixture  may  be  expressed 

e  =  e{v,s,K). 
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Its  differential  is  given  by 


de  =  Tds  -  pdv  +  iidk, 

and  the  condition  defining  the  mixture  adiabat  is 

de  +  pdv  =  Tds  +  p  dA  =  0 . 

Thus,  changes  in  entropy  along  the  adiabat  are  associated  with  chemical  reaction.  In  the  absence  of 
reaction,  the  adiabat  becomes  an  isentrope. 

5.8  Hugoniot.  We  cannot  use  the  mixture  equation  of  state  to  determine  a  Hugoniot  for  the 
mixture  because  the  reactant  isentrope  reference  state  changes  across  a  shock  wave.  In  this  case,  we 
make  use  of  known  equations  of  state  and  Hugoniots  for  the  phases: 

p=  p^{y^,e+Lhg) , 

p  =  Pp(Vp,ep), 

P-Php(V>Pu^^pu)‘ 

And  the  saturation  conditions  still  apply: 

v={^-h)v^+\Vp, 

e  =  (1-A)e,+  Aep. 
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This  is  a  system  of  six  equations  in  eight  variables,  and  the  resulting  Hugoniot  may  be  expressed 


P  =  Ph{v,K\Pu’^u)  ^ 

e  =  eH{v,K]p,,vJ, 

Vr= 

e,=  e^{v,K\p^,v^) , 

Vp=VpH{v,K;p,,vJ, 

ep  =  epHiv^'.Pu>^u)- 

Finally,  a  model  for  the  increase  in  mass  fraction  across  the  shock  is  required: 

A  =  Ay  +  AA^(  v,  e) . 

The  assumption  that  no  reaction  occurs  in  a  shock  wave, 

A  =  Ay, 

is  consistent  with  a  finite  reaction  rate, 
of  the  shock  is  given  by 


5.9  Limiting  Solutions.  The  determinant  of  coefficients  for  the  variational  equations  used  to 
solve  the  mixture  equation  of  state  is  the  same  as  that  given  in  Appendix  B  for  computation  of  the 
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derivatives  of  the  mixture  equation  of  state.  In  the  case  of  thermal  isolation,  a  singularity  arises 
associated  with  the  vanishing  of  this  determinant  at  A  =  0.  In  this  limit,  the  reactant  state  is  identical 
to  the  mixture  state.  Special  considerations  are,  then,  necessary  to  determine  the  product  state.  The 
set  of  equations  obtained  by  differentiating  the  saturation  and  thermal  isolation  conditions  while 
retaining  the  mechanical  equihbrium  condition  is 


(1-A) 


(1-A) 


+  A 

_ £ 

dK 

dK 

\ 

v,e 

\  / 

( 

( ae.'l 

+  A 

p 

v,e 

[aA  J 

v,e 


=  ^r-Sp. 


and 


PsriV>Vri<>erk^^>^D)  -  Ppi^p  .  ©p)  =  0  , 


PsriV’^rl<>ertc*^h^) 


^  dVi '' 


v,e 


( de^ 


=  0 


v,e 


For  vanishing  A,  these  become 


.  ^  , 

\  y  x=o 


=  ^-^PO’ 


^  de,.^ 

.  dK  . 

V  y  x=o 


and 


PpiVpo^epo)  =  Psriv,v^,e^^Ah^) , 


Psriv,v^,ep,^Ah^) 


'dv,'' 

4. 

'ae,' 

A=0 

[ax  J 

=  0 


x=o 
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This  represents  a  system  of  equations  in  the  product  variables  and  the  first  derivatives  of  the  reactant 
variables  with  respect  to  K.  Substituting  the  expressions  for  these  derivatives  into  the  thermal 
isolation  condition  and  coupling  with  the  mechanical  equilibrium  condition  gives  the  equations 


Pp(''pO«®po)  = 
and 


Psri'^’  ^rk  >  “  ^po)  =  ^  ’ 


which  can  be  solved  for  the  product  state  when  X  =  0. 

In  order  to  obtain  a  solution  valid  near  A  =  0,  derivatives  with  respect  to  A  may  be  used  in 
conjunction  with  Taylor's  series.  To  obtain  first  derivatives  of  the  product  variables  and  second 
derivatives  of  the  reactant  variables,  the  set  of  equations  is  differentiated  again,  noting  that 


<^Psr  ;2 


dv 


This  gives 


(1-A) 


(1-A) 


I  aV,  ^ 


/  v,e 


dK^ 

J  v^.e 


I  aA2 


=  2 


/  v,e 


( 

5A 


'  ax  , 

v,e  V  /  v.e 


dh^  j 


=  2 


v,e 


dK 


/  v,e 


dK 


v.e 


^  J 


dK 


v,e 


dK 


v,e 


v,e 


24 


and 


(  \ 
O  Br 

o 

PsriV^^rio^rk^^^D) 

\  5A2  ^ 

+ 

v,e 

r 

<  j 

v,e 

1-) 

v,e 


For  vanishing  A,  these  become 


-  2 

(i'- V  ■ 

<  5A2; 

x=o 

dK 

\ 

{ 

=  2 

(«-«po)  - 

_ p 

k  dk^  J 

A=0 

ax 

\  / 

/  x=o 


+ 

f^pp] 

/ 

CD 

> 

x=o 

CD 

dA 


il{v,v^,e^+Ahij){v-Vpo)  ^  p^iv)  , 


/x=o 


and 


Psriv\v^,e^^Ah[,) 


,  dh^ 


/  A=0 


/  x=o 


-Po{v){v-Vpo)  . 


Substituting  the  expressions  for  the  second  derivatives  into  the  differentiated  thermal  isolation 
condition  and  coupling  with  the  differentiated  mechanical  equilibrium  condition  gives  the  equations 


and 


f^^pl 

+ 

f^ppl 

i^'^pj 

[ax  J 

x=o 

Vn 

[aA  J 

Po('^)  > 


dk 


/  x=o 


1^) 

dh 


/  x=o 


{v-Vpo)[Ps^{v;v^,e^^Lh^)  -|Po(»')]  +  (e- ©po)  ^  eo{v,e) 
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Solutions  for  the  first  derivatives  of  the  product  variables  are,  then, 


eo{v,e)\ 


Sp,] 


/A.0 


de, 


p ) 


-  Poi^) 


^p('^pO*®po) 


and 


A=o 


Poi^)  Psri^’  -  eo{v,e) 


dPp 

K^^PJ 


^p^'^pO>®po) 


The  leading  terms  of  the  Taylor’s  series  expansions  for  the  reactant  and  product  variables  are 

A2  +  .  .  •  , 


X2 


[ 

dv^] 

v^  =  V  +  {V- Vpo)K  + 

I 

1 

P 

[a  J 

A=0 

e,  =  e  +  (e-epo)A  + 


'^P  =  ^PO  " 


'pOf 


V  /  A=0 


l^^p] 

dK 

\ 

x=o 

A  + 


and 


®p  = * 


A=0 


A  + 
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6.  Method  of  Characteristics 


Frontier  uses  characteristic  compatibility  conditions  to  couple  the  Riemann  solutions  for  front 
propagation  to  the  interior  region  solution.  These  must  be  modified  to  account  for  the  effect  of 
reaction  (Grove  1998).  Using  a  local  coordinate  system  (X)  directed  normal  to  the  front,  the 
conditions  applicable  along  streamlines  are 

dt  dX  [  dt  dX  ) 


and 


dt  dx 


A(A,P, 


). 


Here,  h  is  the  enthalpy  defined  as 


h  =  e  +  pv. 


Further,  the  conditions 


dt 


iu±c)^  +  pc 
dx 


i^.{0±c)^ 

dt  dx 


=  [(pc^  -  r p){Vp-  V,)  -  r(ep-  e,)]pA(A,p, 


) 


are  applicable  along  paths  for  which 
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dx 

dt 


=  u±c. 


These  equations  contain  a  source  term  due  to  the  presence  of  chemical  reaction.  In  the  absence  of 
reaction,  the  compatibility  equations  of  gas  dynamics  are  recovered. 

7.  Summary  and  Conclusions 


While  there  is  still  considerable  controversy  surrounding  the  validity  of  the  assumptions  of 
mechanical  equilibrium  and  thermal  isolation,  they  appear  to  represent  the  conditions  prevalent 
during  shock  initiation  of  solid  explosives  with  sufficient  accuracy  for  use  with  a  calibrated  initiation 
model.  Based  on  these  assumptions,  we  have  developed  a  discrete  mixture  equation  of  state  for  use 
with  the  gas  dynamics  conservation  equations  that  takes  advantage  of  the  unique  features  offered  by 
the  front-tracking  method.  We  have  derived  all  of  the  additional  thermodynamic  expressions 
required  to  implement  the  mixture  model  in  the  front-tracking  context.  In  addition  to  incomplete 
equations  of  state  for  pressure  and  temperature,  these  expressions  include  the  sound  speed  and 
Griineisen  functions,  the  temperature  and  specific  heat,  along  with  adiabats  and  Hugoniots.  We  have 
recently  completed  implementation  of  this  approach  in  FRONTIER. 
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Appendix  A: 
Equations  of  State 
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Intentionally  left  blank. 
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A-1.  Equations  of  State  and  Constraints 


A  complete  equation  of  state  has  the  form 


e=e{v,s) , 


such  that 


and 


p{v,s)  =  - 


<  dv  j 


Tiv,s)^ 


^  ae  ^ 
,  as  j 


where  e  is  specific  internal  energy,  v  is  specific  volume,  S  is  specific  entropy,  p  is  pressure,  and 
Tis  temperature.  The  expressions  for  pressure  and  temperature  each  represent  incomplete  equations 
of  state.  By  solving  the  complete  equation  of  state  for  entropy  and  substituting  into  the  incomplete 
equations  of  state,  the  latter  are  represented  as  functions  of  specific  volume  and  internal  energy.  An 
incomplete  equation  of  state  giving  pressure  in  the  form. 


p  =  p{v,e), 

as  required  by  the  gas  dynamics  conservation  equations  can  be  represented  by  a  surface  in  p-v-e 
space.  The  partial  derivative  of  the  equation-of-state  pressure  with  respect  to  internal  energy  is, 
generally,  a  positive  function  of  both  specific  volume  and  internal  energy  denoted  by 


=  Pp(v,e)>0. 

V  j  V 
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Any  constraint  (e.g.,  constant  temperature,  constant  entropy,  etc.)  on  the  processes  that  the 
material  described  by  that  equation  of  state  may  undergo  can  similarly  be  represented  by  such  a 
surface.  The  intersection  of  the  two  surfaces  defines  a  path  along  which  the  constrained  states  must 
lie.  Along  this  path,  each  of  the  variables  (p,  V,  or  e)  may  be  represented  as  a  function  of  either  of 
the  other  two  variables.  For  example. 


e  =  ep,{v), 


and 

P  =  PRiv), 

where  the  subscript  is  chosen  to  identify  the  particular  constraint. 


It  is  possible  to  develop  an  equation  of  state  if  equations  describing  a  constrained  path  and  a 
transverse  derivative  on  that  path  are  known.  For  example,  if  the  derivative  of  pressure  with  respect 
to  internal  energy, 

1^]  = 

\de  )  p 


is  known  on  a  path  as  a  function  of  specific  volume,  an  equation  of  state  valid  near  the  path  is  given 
by 


p{v,e)  =  Pp{v)  +  pj-[v,ep{v)][e  -  ep{v)] . 


The  validity  of  this  expression  is  extended  to  points  arbitrarily  far  from  the  path  if  the  derivative  is 
everywhere  independent  of  internal  energy.  That  is,  the  derivative  must  be  a  function  of  the  specific 
volume  only.  In  this  case,  the  equation  of  state  may  be  written 
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p{v,e)=Pf^{v)  +  p^{v)[e-ep{v)]. 


Using  the  Griineisen  function,  defined  as 


r{v)  =  V 


the  Mie-Griineisen  equation  of  state  is  written 


p{v,  e)  =  Pp{v)  +  [e  -  . 

With  appropriate  reference  and  Griineisen  functions,  it  represents  a  wide  variety  of  commonly  used 
equations  of  state  that  are  linear  in  the  specific  internal  energy.  (A  notable  exception,  offered  by 
Davis  [1993] ,  is  an  equation  of  state  for  detonation  products  that  is  quadratic  in  the  internal  energy.) 

It  is  often  convenient  to  write  these  linear  equations  of  state  in  the  form 

p{v,e)=Pj.{v)  +  pj.{v)e, 


where 


PrC'^)  = 


nv) 


is  a  representation  of  the  Griineisen  function  and 

Pr('')  =  Pni^)  -  Pr('') 


combines  the  reference  and  Griineisen  functions. 
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Solving  for  the  specific  internal  energy  gives 


e{p,v)  =  Bffiv)  + 


y[p-  Pff(V^)] 

r(v) 


P  -  Pr('^) 
Pr('') 


pvj.{v)  -  er(v) , 


where 


Vj.{v)  = 


_J _ 

PrCi^) 


and 


e^{v)  =  Priv)  \^r(v). 


A-2.  Isentropes 

When  the  constraint  is  the  condition  of  constant  entropy,  the  path  is  known  as  an  isentrope.  If 
(Vif,  Gif)  is  a  known  state  on  an  isentrope,  then  the  equation 


e  ^  e[v,${v,^,e,^)]  =  es{v',v,^,e,^) 


describes  the  variation  in  specific  internal  energy  with  specific  volume  along  the  isentrope.  Because 
of  the  fundamental  thermodynamic  identity. 


de=  Tds-pdv, 


the  pressure  and  internal  energy  have  a  special  relation  along  an  isentrope  (ote— 0)t 
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dBs 

dv 


=  =  -PW>es{v]v,^,e^)] . 


This  constitutes  a  first-order  ordinary  differential  equation  for  es{v,  ,  Bi^) .  lu  this  most  general 
form,  it  defies  solution. 

As  we  have  noted,  many  equations  of  state  of  practical  interest  are  linear  in  the  internal  energy. 
In  this  case,  the  differential  equation  can  be  written 


des 

dv 


Pri^)  esiv]v,^,e„)  =  , 


and  a  family  of  solutions  can  be  found  using  an  integrating  factor  (depending  only  on  the  form  of 
the  Griineisen  function)  given  by 


fsiv)  =  exp 


f  Pr('^)dv 


V  ‘'o 


where  V^  is  an  arbitrary  reference  volume.  Defining 


9s('^)  -  Ps(^k)  =  -f  fs(^)  Pri^)  dv. 


the  solutions,  (which  are  independent  of  v^)  are 


^siv-,v,^,e„)  = 


9siv)  -  9siVk)  ^  fs(Vk)ek 
fsiv) 


9siv)  * 
fs(v) 


Thus,  each  isentrope  may  be  described  using  the  single  constant, 
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as(''.  e)  ^  fsiv)  e  -  ggiv)  =  fs{v„)  e„  -  gs{v„)  =  as{v„, e^)=  . 


Defining 


hsiv)  = 


9s(^) 
fsi^)  ’ 


we  may  write 


esiv;Vi^,e„)  =  hs(v)  + 


w 

fsiv) 


[ek-hsiVk)]  =  l^siv)^ 


^$k 


fsiv) 


and 


asiv,e)  =  fsW)[e  -  f7s(v^)]  =  -  /7s(v^^)]  =  as{v,^,e^)  =  . 


If  the  form  of  the  isentrope  is  known  such  that  is{v)  along  with  flfs(  v)  or  hs{y)  can  be  identified, 
an  equation  of  state  can  be  determined  by  differentiating.  Thus, 


Pr('')  = 


1  dfs 
fsiv)  dv 


and 


Pr(»^)  = 


-1  dgs 
fsiv)  dv 


dhs 

dv 


Pr('')  f^siv)  • 


The  equation  of  state  thus  obtained  is  independent  of  the  state  ( ,  61^),  (i.e.,  valid  on  any  isentrope). 
Under  these  conditions,  the  reference  curve  provides  all  the  information  necessary  to  construct  the 
equation  of  state,  including  the  Gruneisen  function. 


We  can  determine  the  pressure  on  the  isentrope  by  substituting  into  the  equation  of  state  or  by 
differentiating  ©sC''  I  ®ir)  • 
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or,  alternatively, 


1  [  9siv)  ^  3,,  dfs 

fsM  [  fs(^)  dv 


Psiv>Vk>e,)  = 


where 


Vsiv)  = 


fsi^) 
Pr('')  ’ 


We  shall  also  have  use  for  the  differential  equation  satisfied  by  pg  {v.Vi^,  e/^) .  It  can  be  obtained 
by  differentiating  the  isentrope  pressure  and  substituting  for  the  isentrope  energy  from  the  equation 
of  state,  giving 


The  integrating  factor  associated  with  this  equation  is  Vg(v) ,  which,  again,  depends  only  on  the  form 
of  the  Griineisen  function. 


We  can  also  relate  the  isentrope  to  the  reference  functions.  Substituting  for  Pj-{v)  gives 


9si^)-9s(^k)  =/  fsi^)  Pr('')efl('')  dv-f  fs{v)Pp{v)  dv. 

Vk 


Noting  that 


=  Pr('')  > 
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the  first  term  of  this  equation  can  be  integrated  by  parts.  Thus, 


9siv)  -  9s(''k)  '  OnW)  -  9M  *  4('')  • 


where 


9Ri^)  -  9R{Vk)  =  -  f  fsi^) 


dv. 


This  gives 


.  ,  /  ^  dRiv)  -  9r{Vi,)  +  fsi^k)[er  ^Ri^k)]  _  .  ^  9r{^ 


9Riv)  + 


for  the  isentrope  energy  and 


,  ,  9r{v)  -  9r(v,)  +  fsiVk)[ek-  eM]  _  „  ^  9^ 

Psi^’^k’^k)  =  Pr^^)  ^ - “ 


9Riy)  + 


for  the  isentrope  pressure. 


The  derivative  of  the  isentrope  pressure  is  the  square  of  the  acoustic  impedance: 


dPs 

dv 


=  Pri^) 


Psi^'^^k’^k) 


It  must  remain  positive. 
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A-3.  Sound-Speed  Functions 


By  differentiating  the  equation  of  state  in  the  form 


p{v,s)=Pj.{v)  +  pr(»^)e(v,s), 


the  square  of  the  acoustic  impedance  can  be  expressed  as  a  function  of  state: 


i^{v,s)  =  - 


'ip]  = 

'  Pj-iv)  _ 

+ 

Pr(i^)  _  dpA 

,  Vr(\/)  dv  ^ 

^  Vr(v^)  dv 

©(v'.s) 


The  square  of  the  sound  speed  and  the  adiabatic  exponent  are  given  by 


c^{v,s)  =  v^i^{v,s) , 


and 


p{v,s) 


These  may  also  be  expressed  as  functions  of  v  and  e: 


i^{v,e)  = 


'  p^{v)  _  dpj.^ 

'  Pr(\^)  _  Qfpr^ 

^  v^{v)  dv  ^ 

^  Vr(\^)  dv  ^ 

e, 


c^{v,e)  =  v^i^{v,e) , 


and 
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vi^{v,e) 

p{v,e) 


These  functions  must  remain  positive  in  order  to  prevent  imaginary  sound  speeds. 


A-4.  Completing  the  Equation  of  State 


A  complete  equation  of  state  can  be  determined  if,  in  addition  to  the  incomplete  equation  of  state 
for  pressure,  an  incomplete  equation  of  state  for  temperature, 

T=T{v,e), 

is  developed  using  suitable  assumptions.  Following  Grove  (1997),  we  solve  the  pressme  equation 
of  state  for  energy,  giving 

e  =  ^  =  pv^iv)  -  e^iv) , 

Pr('') 

so  that  the  specific  internal  energy  differential  is 

dVf  ,  der-  , 

de  =  vAv)dp  +  p—^dv  — ; — dv. 

'  dv  dv 

Substituting  into  the  differential  of  the  complete  equation  of  state, 

de=  Tds-pdv, 


and  rearranging  gives 


ci  der-  . 

Tds  =  vAv)dp  +  p-^[v+  Vj.{v)]dv  -  —^dv. 
'  dv  dv 


Introducing  a  fictitious  pressure,  p^{v),  we  may  write 


42 


Tds  = 


dp„ 


der 


-  -^dv 

dv  dv  dv 


dv 


+  [p-pSv)]—[v*v^{v)]dv+  Vy{y)  d{p - pSv)] 
dv 


The  leading  term  on  the  right-hand  side  of  this  expression  has  been  shown  to  vanish  for 

P»  =Psi^>'^o.>eJf 

where  (V'. ,  e„),  is  a  suitable  reference  state.  Denoting  p-pjy),  referred  to  as  the  stiff  pressure,  by 
p*,  we  may  write 

T{p* ,v) ds  =  p*—[v  +  dv+  vAv)  dp* . 
dv 


Assuming  temperature  in  the  form, 


T{p*,v)  =  T,{p*)Tp,{v), 


gives 


ds 


Vj.{v)]dv 

dv _ 

T,{p*)  Tp.{v) 


+ 


Vr^v)  dp* 
Lip*)  L^iv)  ■ 


Thus, 


P 


■-flv*v,(v)] 

dv 


LiPlL^iv) 


and 
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vAv) 

TAPI  T^-(v)  ' 


Cross  differentiating  gives 


d^S  _  1  dr  /  Xl  d 

P*  I 

1  d 

yrM  ] 

- - - 

dp-dv  T^.iy)  dv  dp- 

TAP-)  dv 

T,Av)  J 

or 


TJpl 


dp' 


[  TAPl 


=  ^  -  w , 


dv 


[v+  \^r('')] 


where  w  is  a  constant,  and  two  differential  equations  result  The  equation  for  T^{p*)  is 


dp* 


and  a  solution  is 


Up*)  = 


p* 

RT^ 


w 

1 


where  RT^  is  a  constant  The  equation  for  7'p.(v')  is 


dv 


(1- w)pr(v^) 


w  dPr^ 
Pfiv)  dv  ^ 


Tp.{v)=0, 


and  a  solution  is 
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^  _  ^oToUsivr-'  [>^s(^)r 


[Pr('^)]‘ 


0  0 


fsi^) 


where  is  a  constant.  Thus, 


T{p,v)= 


fsiv) 


RT. 


W 


and,  substituting  for  pressure. 


T{v,e)  = 


^qTq 

fsiv) 


fsi^) 


(  fsMlPriv)  +Pr{v)  e-Psiv,  VL.  6,0)]  ] 
RT,  p,{v) 


(  [e-esiv]  v^,ej]  /g(y) 
RT^ 


Here,  the  reference  state  (</„ ,  e„),  along  with  the  exponent  w,  may  be  chosen  for  consistency  with 
known  or  assumed  forms  of  the  temperature.  Solving  for  internal  energy  gives 

RT^  (  fsiv)  T  Y''^ 

©(T’.v')  = -rrrl I  +es{v;v^,ej, 


fsiv) 


V  '■o  T^j 


and  the  specific  heat  at  constant  volume  is 

Cy{T,v)  = 


'  ae ' 

R 

fsM  T  1 

,aT, 

V 

r,  7-J 

This  is  inconsistent  with  results  from  quantum  statistical  mechanics,  which  show  that  Cy  has  no 
volume  dependence  and  is  bounded  for  large  temperatures  ( Vincenti  and  Kruger  1967;  Callen  1 960) . 
However,  for  polyatomic  gases  at  low  and  very  high  temperatures  and  for  solids  at  very  high 
temperatures,  Cy  is  nearly  constant,  corresponding  to  =  1 .  In  this  case. 
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and  the  temperature  is  given  by 


Pci^)  = 

or 

ec{v)  =  esiv;v„,ej, 

where  i^)  specific  interatomic  potential  energy.  As  a  practical  alternative,  may  be 
determined  with  reference  to  a  known  temperature.  If  Tj  is  the  known  temperature  at  the  state 
(v^,  Bj) ,  then  the  expression 

3^=^srW^vo'^j 
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constrains  the  state  ,  e„ ) .  In  this  case,  the  proper  behavior  at  cold  temperatures  may  not  be 
realized. 

The  temperamre  on  an  isentrope  is  given  by 


c^oPr(^) 


'VO 


Specific  entropy  in  the  form. 


/ 

s{p\v)  =  c Jog 

\ 


P*Vsi^) 


\ 


) 


I 


satisfies 


P*  —  [v+Vj.{v)] 

dv _ 

T{p\v) 


Oyp 

Vs(v)  dv 


and 


as  \  ^  M'')  ,  £» 
dp- 1,  T(p-,v)  p- 


but  violates  the  Nemst  postulate.  The  complete  equation  of  state  is 


e(v,s) 


fh. 

fs(^) 


+  Bciv) • 
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A-5.  Hugoniots 


We  have  shown  that  an  isentrope  can  be  used  to  determine  an  equation  of  state.  A  shock 
Hugoniot  is  another  path  in  the  p-\^  plane  that  is  used  in  equations  of  state.  Rather  than  describing 
a  process,  it  represents  the  locus  of  all  possible  states  on  the  downstream  side  of  a  shock  wave  for 
a  specified  upstream  state. 

If  U  denotes  the  shock  propagation  velocity  and  u  the  component  of  the  particle  velocity  normal 
to  the  shock  wave,  the  equations  expressing  conservation  of  mass,  momentum,  and  energy  across 
the  shock  wave  may  be  written 

^^v{U-uJ  =  v^A"u, 

{U-uJA^u  =  A^pv,, 


and 

A“e  +  A“(p  v)  =  {U-  u^)  A“t/  -  . 

Here,  unsubscripted  variables  refer  to  the  downstream  state,  and  the  subscript  u  denotes  known 
conditions  on  the  upstream  side  of  the  shock  wave.  The  operator  A"  represents  changes  across  the 
shock: 

A"Q  =  £7-qry. 

For  convenience,  we  also  define 

"Av'=  v^-v. 


Eliminating  U  -U^  between  the  mass  and  momentum  conservation  equations  gives 


(A^uf  =  A‘^p  “Av. 
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Substituting  the  expressions  for  {U  -  u^)  L"u  and  (A"t;  Y  into  the  energy  equation  gives  the 
Hugoniot  relation: 

A“e  =  (p„+|A“p)  “Av. 

This  represents  a  surface  in  p-v-e  space.  The  intersection  of  this  surface  with  the  equation-of-state 
surface  is  obtained  by  inverting  the  equation  of  state  to  give  specific  internal  energy  as  a  function 
of  pressure  and  specific  volume  and  substituting  into  the  foregoing  equation.  Noting  that  the 
upstream  material  may  be  described  by  a  different  equation  of  state  (e.g.,  due  to  chemical  reaction), 
this  gives 

e{p,v)  -  ejp^,vj  =  {p^+  ^A'^p)  “Av, 

which  represents  a  Hugoniot  in  the  p-»^  plane.  In  principle,  this  can  be  written 

P  =  PHi^’Pu’^u)^ 


or,  alternatively, 


A“p  =  A“p„(\/;p„,vJ. 


The  Hugoniot  energy  is  related  to  the  Hugoniot  pressure  by  the  Hugoniot  relation; 


eni^'^Pu’^u)  =  e^iPu>^u)  ^  j[Pu^ Pni^^Pu’^u)'^ 


or 


N^e^iv; p^,v^)  =  [p„+  ^A^Pniv, P„,  vj]  "Av^. 


This  may  be  differentiated  to  give 
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±£!  =  if  -  p„-  p„{v.p,.vj  I  =  if  -  A"p„(»';p„.i/„)  -  P„  • 

dv  2  dv  \  ^\  dv  j 


For  an  equation  of  state  that  is  linear  in  internal  energy,  we  may  write 


Ph(^’Pu>^u)  = 


Priv)  ^  n  -  Vu)]Pu  ^  Pr('^)  ^u^Pu  >  ^u) 


A“p„(v^;p„,\/u) 


Pr('^)  -  P»-^Pr(^)tgt/(P»y»0  ^ 


for  the  Hugoniot  pressure.  Here  the  Hugoniot  factor  is  defined  as 

Hugoniots  are  valid  in  compression  (V'  ^  vj  only  and  may  be  thought  of  as  emanating  firom  the 
upstream  state.  If  the  Hugoniot  factor  vanishes,  there  may  also  be  a  lower  limit  on  the  specific 
volume  that  may  be  obtained  by  shock  compression.  Thus, 


V  “Av 


Segletes  (1991)  has  identified  this  inequality  as  a  stability  requirement  for  the  Hugoniot-referenced 
Mie-Gruneisen  equation  of  state.  The  minimum  specific  volume  is  obtained  by  solving  the  equation 


Pr^^min ) 


V  -V  ■ 
min 
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The  Hugoniot  pressure  and  energy  become  unbounded  as  the  specific  volume  approaches  this 
minimum.  This  singularity  is  due  only  to  the  form  of  the  Griineisen  function. 

The  Hugoniot  and  isentrope  functions  for  an  equation  of  state  that  is  linear  in  internal  energy  are 
related  by  the  expression 

Pu>'^u)  ^  Pu^^u)]  ■ 

Differentiating  gives 


^  Pfi^)  ’  Pu^^u)  -  Psi^’^k’^k)  ^  Ph(^’  Pu’^u)^ 

av  av  2 

dPr 

av 


Thus,  the  derivative  of  the  isentrope  at  the  Hugoniot  is 


dPs 

dv 


H 


dP„ 

dv 


+  ^pf{v)A^PH{v,Pu,Vu). 


Evaluated  at  the  upstream  state,  this  becomes 


dPs 

dv 


dPn 

dv 


> 


U 


showing  that  the  Hugoniot  and  isentrope  are  tangent  as  they  pass  through  that  state. 
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A-6.  The  Griineisen  Function 


The  Griineisen  function  is  commonly  represented  in  the  linear  form  (Steinberg,  Cochran,  and 


Guinan  1980) 


so  that 


r(u)  =  r„+Pr/. 


The  integrating  factor  for  the  isentrope  energy  becomes 


fgiv)  =  exp  J Pj-{v)dv  =  (v/v’o)^°exp[prJ\^-Vj,)] 


The  derivative  of  Pr('^)  is 


dPr  ^0 


and  the  square  of  the  acoustic  impedance  may  be  written 


,  f  rj  _  dp,  Pr,  r„(r.*i)  „ 

i^{v,e)  =  Pr  •—  PrW)  -  ^  t  — i  *  — —  - - —  e 

'o  V  dV  Vr  V 


The  complete  equation  of  state  is 


e{v,s)  =  RTJ—  exp 


v„  V  °  s  V 

-2-  exp  — - 
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where 


With  Pro=0,  the  Griineisen  function  is  constant,  and,  with  ro=0,  it  is  proportional  to  the  specific 
volume.  We  shall  find  it  useful  to  specialize  some  of  the  expressions  derived  in  the  foregoing 
sections  to  these  cases. 


A-6.1  Constant  Gruneisen  Function.  The  simplest  form  of  the  Griineisen  function  is  the 
constant-f  form  (P[-o=0): 


Pr('')  = 


V' 


Differentiating  gives 

=  -lo  . 

dv 

The  isentrope  energy  and  pressure  integrating  factors  become 


fs{v)={v/v^f\ 


and 


Vsiv)  = 


fsi^)  ^ 


The  square  of  the  acoustic  impedance  is 


dPr 

dv 


+  — — - - 6 


53 


The  Hugoniot  factor  and  minimum  specific  volume  are 


=  1  - 


> 


and 


V  ■  =  — - — 
mm  p  ^  2 


The  complete  equation  of  state  is 


e{v,s)  =  RT^ 


To 

exp 

+  es{v]v^,ej. 


A-6.2  Proportional  Griineisen  Function.  Another  simple  form  of  the  Gruneisen  function  is 
the  proportional  form  (Fq  =  0): 


Pr('^)  = 


IM 

V 


where  Pro  is  a  constant  and  the  derivative  of  Pr  (v')  vanishes.  The  isentrope  energy  and  pressure 
integrating  factors  become 

/s(v)  =exp[Pr^(v-\^o)]> 


and 

(V)  =  .  exp[Pr^(y-yJl 

®  PrC'')  Pr, 

The  square  of  the  acoustic  impedance  is 
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i^{v,e) 


|^Pr,Pr('')- 


dv 

/ 


The  Hugoniot  factor  and  minimum  specific  volume  are 


and 


min 


The  complete  equation  of  state  is 


/ 

e{v,s)  =  RT^exp 

< 


^es{v;v„,ej. 


A-7.  Polytropic  Gas  Equation  of  State 


A-7.1  Incomplete  Equation  of  State  for  Pressure.  With  a  constant  Griineisen  coefficient, 
f"o  =  Yo  reference  to  zero  pressure  and  internal  energy,  the  polytropic  gas  equation  of  state 

results: 


p{v,e)=^^^e. 

V 

The  function  Pp  (v)  may  be  identified  as 
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Pr(’^)  = 


Yo-1 

V 


> 


and  Pr(v')  vanishes. 

A-7.2  Isentrope.  Thus,  the  integrating  factor  for  the  isentrope  energy  is 


=  exp 


=  exp[(Yo-1)  log(W\/o)]  =  . 


The  function  Qsiv)  vanishes.  The  internal  energy  and  pressure  on  the  isentrope  are,  then,  given  by 


W) 


and 


W 

fsi^) 


^kPri^)  = 


^l2Zle^{Vi^/vy°  =  p^{v^/v)'^°. 

^k 


Thus,  on  an  isentrope. 


and 


V 


Yo-1 


e  = 


Pv'^°=Pk^l''  ■ 
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A-7.3  Sound-Speed  Functions.  The  square  of  the  acoustic  impedance  is 


/2  = 


^  Pri^) 


Yo(Yo-1)  ^ 
— — - —  e  = 


the  square  of  the  sound  speed  is 

c^=  v^^/^-Yo(Yo-1)e  =  YoP''. 

and  the  adiabatic  exponent  is 


These  are  positive  for  positive  specific  internal  energy. 


A-7.4  Complete  Equation  of  State.  The  temperature  of  a  polytropic  gas  is  known  to  vanish 
for  vanishing  specific  internal  energy.  This  is  consistent  with  a  vanishing  specific  interatomic 
potential  energy; 

eciv)  =  es{v;v„,ej=0, 

which  implies  that 


and 


Thus,  the  stiff  pressure  is  simply  equal  to  the  pressure,  and  the  temperature  is  given  by 

T-^=E^  =  ^-£le. 

Cyo  ^  FI 
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The  integrating  factor  for  the  isentrope  pressure  is 


Vgiv)  = 


^s(^) 

Pr('') 


Yo-1 


The  entropy  may  then  be  written 


f 

S  =  CjOQ 


pyivlvj’''' 

(Yo-drt; 


\ 

=  cjog 

^  V  ^ 

p  V^° 

=  CjOQ 

Yn"1  ^ 

V  °  e 

Yo 

,,Yo-1  . 

/ 

[Po^o  i 

k  ^0  ; 

9 


where 


eo-RT, 


and 


e. 


Po=(Yo-1)  — 


I 


and  the  complete  equation  of  state  is 


Yo-1 

e  =  RT, 

exp 

k  ) 

A-7^  Hugoniot.  The  Hugoniot  factor  is 


Yo-1 

2v 


and  the  Hugoniot  pressure  is  given  by 
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PHi^>Pu>^u) 


(Yo+1)v^»-(Yo-1)v^ 
(Yo+1)''  -(Yo-I)'',;^" 


The  minimum  specific  volume  is 


V  ■  = 
mm 


A-8.  Stiffened  Polytropic  Gas  Equation  of  State 


A-8.1  Incomplete  Equation  of  State  for  Pressure.  The  polytropic  gas  equation  of  state  is 
“stiffened”  by  the  use  of  a  constant  pressure  reference: 


Pr{v)  =  YoP»  , 


so  that 


p(v,e)  =  YoP„ 


e. 


V 


The  function  pr  {v)  remains  unchanged  while 


Pj.{v)  =  . 


A-8.2  Isentrope.  The  integrating  factor  for  the  isentrope  energy  remains  unchanged  and 


9si^)  =  -P-''o(''/''o)^“  • 


The  internal  energy  and  pressure  along  the  isentrope  are,  then,  given  by 
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’  ^k)  =  {yjv)'^°  -pjv, 

and 

Ps(^’^k’^k)  =  P»+(Yo-1)(Pco+e^/''/f)K/'')^®- 


Thus,  on  an  isentrope, 


{p«+  e!  v)  /°  =  (p„+  e^l  v^)  vl° , 


and 

(P-PJ  v''°  =  iPk- pj  vl\ 


A-8.3  Sound-Speed  Functions.  The  square  of  the  acoustic  impedance  is 


PrM  _ 

dpr^ 

+ 

Pr('') 

_  dp,) 

e=Vo(y,-1)r  eUVo-1f 

dv  J 

dv 

> 

V  ''i  V 

the  square  of  the  sound  speed  is 

c2=  v^i^  =  \^{y^-^){p„v+e), 


and  the  adiabatic  exponent  is 


Yo(Yo-1)(P.o<^+ e) 

pv  Yo(P.^+®)  +  e 


60 


A-8.4  Complete  Equation  of  State.  Choosing  (K. ,  e„ )  so  that 


=  -P„ 


associates  with  the  “cold  pressure”: 

Pci^)  =  . 

The  resulting  specific  interatomic  potential  energy  is 


eci^)  =  es{v;v^,ej  =  -p„v . 


Thus,  the  stiff  pressure  is 


P*  =  P-P.  . 


and  the  temperature  may  be  written 


j_  p*v  _  ^  P~^ 


R 


'VO 


The  integrating  factor,  Vq{v),  for  the  isentrope  pressure  is  unchanged  and  the  specific  entropy 
is  given  by 


s{p,v)  =  c Jog 


(  ip-pjv{v/vy°-^] 


(Yo-D^t; 


=  cjog 


'(p-pJv-^°' 
V  {Po-PJvJ°  J 


or 


s(\^,e)  =  c^Jog[ 


u^°'^(e  +p„v 


I 
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where 


j 


and 


Po  =  YoP-  +  (Yo-1) 


> 


and  the  complete  equation  of  state  is 


Yo-I 

e-RT, 

V 

\  J 

exp 

A-8.5  Hugoniot  The  Hugoniot  factor  and  the  minimum  specific  volume  remain  unchanged 
from  the  polytropic  gas  case,  and  the  Hugoniot  pressure  is  given  by 


[(Yo^1)'^a-(Yo-1)*^]Pu-2YoP,K-^) 

(Yo+1)''  -(Yo-1)''«; 


A-9.  JWL  Equation  of  State 

A-9.1  Incomplete  Equation  of  State  for  Pressure.  In  the  case  of  the  JWL  equation  of  state, 
commonly  applied  to  detonation  products,  we  begin  with  a  knowledge  of  the  isentrope  emanating 
from  the  Chapman-Jouget  state,  (Vj,  e,).  Lee  and  Homig  (1969)  chose  the  form 
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es{v]Vj,ej)  =  hgiv)  + 


ej  -  hsivj) 

{vivjr 


so  that 


dhs  ^  0)  ej-hsivj) 
dv  V  {v/Vj)'^ 


where 

hsiv)  =  eiexp(-piV)  +  e2exp(-p2u) , 
and 

dh^ 

=  p^exp{-p^v)  +  p^expi-p^v) . 


e„  .  ^  m=1,2, 

Pm 

and  Pi,  P2,  Pi,  P2,  and  (1)  are  constants  that  are  determined  with  reference  to  cylinder  expansion  test 
pressure- volume  data.  By  comparison  with  our  general  solution  for  the  isentrope  energy,  we  see  that 
this  corresponds  to  a  constant  Gruneisen  coefficient,  =  CO,  and 


Thus, 


Pr('^)  = 


IM  = 

V 


1  dfs 

fsiv)  dv 


V 
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and 


Pr('^)  =  - 


dhs 

dv 


/ 

-  Pr('')  =  Pi 


/ 

exp(-piV^)  +P2 

\ 


exp(-p2v) , 


where 


/77=1,2  . 


The  JWL  equation  of  state  can  be  written 


p{v,e)  =  Pi 


[1-^] 

exp(-pi\^)  +  P2 

''J 

'' J 

exp(-p2v)  +-^e  . 


A-9.2  Isentrope.  The  general  isentrope  has  the  same  form  as  the  Qiapman-Jouget  isentrope: 


es(v,v^,e^)  =  Pg(v^)  + 


e„-hs{v„) 

{viv.r 


or 


Ps(^>^k>ek)  = 


dhs  ^  gj  e„-hs{v^) 
dv  V  [vlVf^^ 


A-9.3  Sound-Speed  Functions.  Noting  that 


exp(-pi\^)  +  P2 


exp(-p2v)  , 


we  can  obtain  an  expression  for  the  square  of  the  sound  speed: 
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c^{y,e)  = 


a)(o)+1) 

ip^vf  _  A 

exp(-piv)  +  ©2 

{P2Vf  _  J 

expC-PgV)  +  e 

a)((o+1)  j 

to  ((0+1)  J 

The  square  of  the  acoustic  impedance  and  the  adiabatic  exponent  may  then  be  determined  from  their 
definitions. 


A-9.4  Complete  Equation  of  State.  The  temperature  may  be  written 


r= 


^  ^vo 


e-es{v;v^,ej 


'VO 


) 


where  and  e„  are  determined  with  reference  to  a  known  temperature.  The  integrating  factor  for 
the  isentrope  pressure  is 


V's('')  = 


fsiv) 

Pr('') 


=  -{v/vj 


0) 


(0 


Thus,  the  specific  entropy  is 


s{p,v)  =  c Jog 


[p-Psiv,v^,eJ]v{v/vJ^ 


0)2  c  7 


VO  o 


=  CjOQ 


'  [p  -Ps(^  ;v;..0] 

,  [Po-PsC'^o^'C-eJ]  ) 


or 


s(v^,e)  =  c„Jog 


'  y^[e  -egiv  ;v^,ej]  ^ 


where 
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e^  =  RT^  +  es{v^;v^,ej 


and 


e„ 


Po  =  Pr(0-"“  — 


and  the  complete  equation  of  state  is 


f 

(D 

e{v,s)  =  RT^ 

V 

\ 

exp 

+  es{v,v^,ej 


A-9.5  Hugoniot  The  Hugoniot  pressure  is  given  by 


PHi^’Pu’^u)  = 


2[Pf-{v) +  pj-{v)eJv  +  (o"Av 

2  V'  -  0)  “Av^ 


9 


and  the  minimum  specific  volume  is 


(0 

1/  .  =  ■  —  V 

^mtn  _  .  o  u 


aj+2 


A-IO.  A  JWL  Alternative 


A-lO.l  Incomplete  Equation  of  State  for  Pressure.  An  equation  of  state  with  a  constant  T  is 
typical  of  gases.  Although  Lee  and  Tarver  (1980)  have  used  this  form  of  the  JWL  equation  of  state 
to  represent  solid  materials,  they  often  exhibit  behavior  consistent  with  a  proportional  Griineisen 
function: 
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Pr(v')  = 


r(v) 

V 


where  is  a  constant.  An  alternative  form  of  the  JWL  equation  of  state  reflecting  this  behavior 
may  be  obtained  from  an  isentrope  in  the  form, 


esiv]v^,e;)  =  hsiv)  + 


er-hs{Vr) 

exp[prji^-\^,)] 


so  that 


Ps{v]v^,e;)  = 


dhs  Prje,- /)s(v^p)] 

dv  exp[ppji^-»/,)] 


where  hs{v)  remains  as  in  the  constant- F  JWL  formulation  and  e^)  represents  the  state  defining 

the  reference  isentrope.  In  this  case, 

fs{v)=exp[pj.J.v-v^)]. 


Thus, 


Pr('')  =  Pr„. 


and 

Pj-iv)  =  ©^(Pi  -  Pp^)exp(-Pi\/)  +  e2(p2-  Pro)®^P(”p2'')  • 


The  equation  of  state  is 

p{v,e)  =  ei(Pi  -  Pr^)exp(-piV)  +  e2(p2-  Pr„)®xp(-p2'')  +  Pr,,®- 
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A-10.2  Isentrope.  The  isentrope  through  the  state  (14,  e*.)  is 


Gk-hsiVk) 

exp[prj\/-\/;,)] 


I 


or 


Psi^’^k>ek)  = 


dhg  ^  Pri^k-  f^si^k)] 

dv  exp[prJ»^-V4)] 


A-10.3  Sound-Speed  Functions.  With 

=  pi(pr  -  pi)exp(-piV)  +P2(Pr -P2)®XP("P2'')> 

dv  ° 

the  square  of  the  acoustic  impedance  is 

i^{v,e)  =  e^{p^-  Pr^)(Pi  +  Pr^)exp(-piV)  +  egCpg-  Pr^)(P2  +  Pro)®^P(~p2'')  Pr^®- 

The  square  of  the  sound  speed  and  the  adiabatic  exponent  may  then  be  determined  from  their 
definitions. 

A-10.4  Complete  Equation  of  State.  The  temperature  may  be  written 

7=  [p-Ps(^:K..g«)3  ^  e-eg(y;t4,e^) 

Pfo  ^vo  ^VO 
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where  v„  and  are  determined  with  reference  to  a  known  temperature.  For  a  proportional 
Griineisen  function,  the  isentrope  pressure  integrating  factor  is 


exp[pr 

_  0 _ 

Pr 

'  0 


and  the  specific  entropy  is 


/ 

S{p,v)=C  Jog 


[p  -  PsiV’  eJ] exp[pro(y- Vq)] 


J 


( 


[p  -Psiv  ;y.,eJ]exp(Pr/) 
IPo  -  Ps('^o ;  K. .  e„)  ]  exp  ( pj-v^) 


\ 


J 


or 

[e  -esiv  ;\/,,eJ]exp(Pr/)  ^ 
[@0- es(''o;''-.eJlexp(Pr^vJ  ^ 

where 

=  es{v^;v^,ej 
and 

Po  =  PrK)^Pr„eo> 


s(v^,e)  =  c^Jog 


and  the  complete  equation  of  state  is 


/ 

e{v,s)  =  RT^exp 


+  es{v;v„,ej 


'0  / 


The  Hugoniot  pressure  is  given  by 
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9 


L'^Pniv)  = 


and  the  minimum  specific  volume  is 


VL  i„=  v-2Vr 
min  u  I  o 


Measurement  of  a  solid  isentrope  for  the  purpose  of  calibrating  these  equations,  however, 
presents  some  difficulty.  The  equation  is  easier  to  calibrate  with  reference  to  Hugoniot  data.  This 
is  facilitated  by  converting  the  measured  shock  and  particle  velocities  to  pressures  and  specific 
volumes  using  the  expressions 


and 


^Av  = 


The  Hugoniot  expression  may  then  be  fit  to  this  data  in  order  to  obtain  values  for  p, ,  P2,  Pi ,  and  P2. 
The  Griineisen  constant,  Pr,  may  also  be  included  in  the  fit  or  may  be  considered  known  from  other 
sources.  In  general,  this  procedure  does  not  provide  good  fits  to  linear  Hugoniots. 


A-11.  Hugoniot  Equations  of  State 


A-ll.l  Hugoniot  Reference  Functions.  The  ideal  gas  and  JWL  equations  of  state  have  been 
shown  to  have  the  general  form  of  the  Mie-Griineisen  equation  of  state.  However,  the  latter  is  most 
closely  associated  with  the  use  of  a  shock  Hugoniot  as  a  reference.  An  experimentally  determmed 
Hugoniot  can  be  used  as  the  basis  for  developing  an  equation  of  state.  Generally,  however,  p-v 
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measurements  are  not  made.  Rather,  shock  and  particle  velocities  are  measured  in  impact 
experiments  and  a  curve  fit  in  the  form 

is  generated  for  use  in  lieu  of  the  energy  equation.  Here,  ambient  conditions,  denoted  by  the 
subscript  a,  give  the  upstream  state.  The  limit  of  the  shock  speed  as  the  upstream  state  is  approached 
is  equal  to  the  ambient  sound  speed,  Cg.  Thus, 

U„{0)=Cg. 


In  order  to  be  of  practical  use,  the  Hugoniot  must  be  expressed  in  terms  of  pressure  and  specific 
volume,  making  use  of  conservation  of  mass  and  momentum.  Substituting  into  the  momentum 
equation  gives 


^Av 


The  functional  form  of  UniA^u)  must  be  chosen  such  that  this  equation  can  be  solved  for 
A°Ph{v,  Pg ,  Vg).  A  problem  arises  in  generalizing  the  equation  of  state  to  arbitrary  upstream  states 
when  the  experimental  Hugoniot  is  used  since  it  is  not  clear  to  what  extent  L/^(A^t/)  depends  on  the 
upstream  state  beyond  its  known  hmiting  value. 


For  example,  the  simplest  form  of  the  reference  Hugoniot  is  associated  with  a  linear  fit  of  shock 
velocity  as  a  function  of  particle  velocity  using  data  for  material  shocked  from  the  ambient  state: 

U=Cg+SgA^U. 


Many  solid  materials  conform  closely  to  this  behavior.  The  consequent  Hugoniot  pressure  and 
energy  jumps  are 


A^PHW;Pg,Vg)  = 


Cg^AV 

{Vg-Sg^AVf 
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and 


A^eH{v]p^,v^)  = 


Pa^ 


cf  ^Av  ^ 


^Av  . 


2{v^-S^^Avf) 

Note  that  the  reference  Hugoniot  (and  the  consequent  equation  of  state)  are  singular  at 

S.-1 


MIN 


V^. 


This  singularity  arises  out  of  the  reference  function  choice  and  is  distinct  from  the  singularity  in  the 
general  Hugoniot  associated  with  the  Gruneisen  function,  which  occurs  at  v  =  It  limits  the 
utility  of  the  linear  Hugoniot  equation  of  state  to  v>  Vmin,  even  though  the  material  it  describes  can 
be  compressed  to  smaller  volumes  by  isentropic  processes. 

The  isentrope  and  sound-speed  function  expressions  make  reference  to  the  derivatives  of  the 
Hugoniot  functions.  These  are 

dPH  ^  ^  Va*s^^Av  A^Ph  ^ 

dv  Vg- Sg^Av  ^Av 


and 

dv  {Vg  -  Sg  ^Av)^  V  ^  ^a~ 


A-11.2  Incomplete  Equation  of  State  for  Pressure.  The  equation  of  state  with  reference  to 
the  ambient  Hugoniot  can  be  written 

p{v,e]Pg,Vg)  =  PH('';Pa.'^a)  ■"■^~t®“®H^'^>Pa’''a)]  « 


or,  substituting  for  the  Hugoniot  energy  from  the  Hugoniot  relation. 
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p{v,e] p^,v^)  =  f^iv-.v^) Puiv'. p^,v^)  +  ^  {e -  e^- . 

In  writing  these  equations,  we  have  emphasized  the  dependence  on  the  choice  of  an  upstream  state. 
The  ambient  state  is  commonly  used  and  is  a  good  choice  for  describing  a  material  that  has  been 
shocked  from  that  state.  More  generally,  materials  may  be  multiply  shocked  and  the  Hugoniot 
reference  for  the  most  recent  upstream  state  should  be  more  accurate.  However,  because  of  the  way 
Hugoniots  are  measured,  such  generality  is  seldom  achieved. 

The  function,  Pf{v),  is  given  by 

Pj-{v] p^ , vj  =  Ph{v;p^ ,v^)-  pr(v)  p^,v^), 


or,  eliminating  the  Hugoniot  energy. 


PAv.P^  .  ''a)  =  fni'' Pd''’  Pa  ■  ''a)  -  Prl'')  (®a 


Differentiating  gives 


dv 


dPn  ,  ,  deu  ,  ,  dpj- 

-T  ■  Pr 0^)  -r  ~  ’ 

dv  dv  dv 


or 


dPr 

dv 


[P3+  Ph{v;  Pa,  uJ]  -  e^iv]  p^,v^) 


Gfpr 

dv 


A-113  Isentrope.  Recalling  the  relationship  between  the  isentrope  and  reference  energy 
functions,  and  using  the  differentiated  Hugoniot  relation,  we  may  write 
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9Hi^’  Pa^^a)-  dni'^k  >Pa>^a)*  ^s(^/r)  [g/r"  ^  Pa  >  ^a  )3 

W)  ’ 


where 


Pa’  ^a)  -  Qn^^k’Pa’  '^a)  =  "  ^  /  W 


A^Pni^’Pa’^a)  + 


dPn 

dv 


'Av 


dv. 


In  general,  the  integral  must  be  evaluated  numerically.  Forest  ( 1 995)  has  given  a  method  to  facilitate 
the  integration.  A  principal  disadvantage  of  the  linear  Hugoniot  is  the  absence  of  analytic  expressions 
for  the  isentrope  internal  energy  and  pressure. 

A-11.4  Sound-Speed  Functions  and  Complete  Equation  of  State.  The  sound  speed  functions, 
temperature,  and  complete  equation  of  state  are  determined  by  substitution  into  the  appropriate 
general  equations.  No  special  insight  is  afforded  by  this  exercise. 

A-IU  Hugoniot  The  pressure  on  a  Hugoniot  with  arbitrary  upstream  conditions,  denoted  by 
the  subscript  U,  is  given  by 

PH('^>Pa>^a>Pu’^u)  = 

- - {fH{v;v^)  Ph{V]  p^,v^)  -  [\-  fH{v-,v^)]p^^[\-  fH{V,V^)]Pu^  Pr(^)  i^u-  • 
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This  Hugoniot  is  valid  for 


^MIN^  '^min  ’ 

where  depends  on  the  form  of  the  Gruneisen  function. 
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Intentionally  left  blank. 
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Appendix  B: 

Mixture  Equations  and  Derivatives 
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Intentionally  left  blank. 
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B-1.  Equations 


The  equations  required  to  treat  mixtures  in  mechanical  and  thermal  equilibrium  are  as  follows 
Conservation  Equations: 

+V-(pt/)  =0, 

dt 

—  (pti)  +V-(pat/)  +Vp  =  0, 

dt 

■|-(pe7.)  +  V-(pe7-t/)  +  V-(pw)  =  0, 

e  =  e  +  -U-U, 

'  2 

—  (pA)  +  V'(pAu)  =  pA(A,p . ) . 

dt 

Mixture  Equations: 

v=i^  -K)v^+KVp, 
e  =  (1-A)e,+  Aep, 
p=  p^{v^,e+ Ah  d), 

P  =  Pp('^p>ep)> 
r=  T^{v^,e^+Ahg) , 

Tp{Vp,ep). 
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The  equations  required  to  treat  mixtures  in  mechanical  equilibrium  and  thermal  isolation  are  as 
follows: 

Conservation  Equations: 

iP-  +  v-(pt/)  =  o, 

dt 

—  (po)  +v-(puu)  +vp  =  0, 
dt 

—  ipej)  +  V-{pejU)  +v-(pt/)  =0, 

dt 

eT=  e  +  -u-u, 

—  (pA)+V-(pA£i)  =  pA(A,p . ), 

dt 

A(p^j  +  V(pi^,^u)=0. 
dt 

^(perfc)+V-(pe^u)  =  0. 

dt 

Mixture  Equations: 

v={^  -K)v^+KVp, 

e  =  (1 -A)e^+ Ae^, 

P= 

p  =  Pp{Vp,ep), 

P  =  Psr(V’Vri,’^rk^^^D)- 
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In  neither  case  is  closed-form  solution  of  the  mixture  equations  possible.  It  is  necessary  to  solve 
the  equations  using  an  iterative  approach  (e.g.,  Newton's  method). 

B-2.  Derivatives 


ha  order  to  determine  many  thermodynamic  functions  of  the  mixture,  derivatives  of  the 
dependent  state  variables  with  respect  to  v,  e,  and  A  are  required.  A  general  expression  for  the 
differentiated  mixture  equations  that  applies  to  both  thermal  equilibrium  and  thermal  isolation  is 
given  by 


1-A 

0 

A 

o 

1  . 

cT 

0 

1-A 

0 

A 

^2 

C2 

^31 

^2 

®33 

^34 

^3 

0 

,^41 

^42 

^43 

^44 

0 

Here,  the  a^j  and  a^y  coefficients  differ  depending  on  the  model.  For  thermal  equilibrium,  they  are 


^i- 


'ap/ 

l^Pp] 

[aej 

ia®pj 

''p  . 

and 


El 


E 
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For  thermal  isolation,  they  are 


0 


(  ^Pp 


V 


and 


1  0  o]-[p,  1  0  o].[Pp  1  0  0], 


where  the  square  of  the  acoustic  impedance  is 


dv. 


Sr 


The  derivative  vector,  bf ,  and  the  constant  vector,  C, ,  differ  depending  on  whether  derivatives  are 
taken  with  respect  to  v,  e,  or  K.  For  derivatives  with  respect  to  v,  they  are 
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For  derivatives  with  respect  to  6,  they  are 
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For  derivatives  with  respect  to  A,  they  are 
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and 

(^i=[^r-^p  0  0]. 

The  determinant  of  coefficients,  then,  is 
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where 
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Solving  for  derivatives  with- respect  to  v/ gives 
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and 


Similarly,  solving  for  derivatives  with  respect  to  e  gives 
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Derivatives  with  respect  to  A  are  given  by 
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Considerable  simplification  occurs  in  the  case  of  thermal  isolation,  where 
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The  determinant  of  coefficients  becomes 
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Derivatives  with  respect  to  v,  then,  are 
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Derivatives  with  respect  to  e  are  given  by 
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